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We explore a possible phase where chiral SU{2)l x SU{2)r symmetry is spontaneously broken 
while its center Z2 symmetry remains unbroken and its consequence on thermal quantities. In this 
phase, chiral symmetry breaking is driven by a quartic quark condensate although a bilinear quark 
condensate vanishes. A Ginzburg-Landau free energy leads to a new tricritical point (TCP) between 
the Z2 broken and unbroken phases. Furthermore, a critical point can appear even in the chiral limit 
where explicit breaking is turned off, instead of a TCP at which restoration of chiral and its center 
symmetries takes place simultaneously. The net quark number density exhibits an abrupt change 
near the restoration of the center symmetry rather than that of the chiral symmetry. Hadron masses 
in possible phases are also studied in a linear sigma model. We show that, in the Z2 symmetric 
phase, the qq-type scalar meson with zero isospin I — splits from the gg-type pseudoscalar meson 
with I — 1. 

PACS numbers: ll.30.Rd, ll.30.Ly, 25.75.Nq, 21.65.Qr 
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1. INTRODUCTION 

Properties of hot and/or dense QCD matter has been 
extensively studied within chiral approaches [l|. Our 
knowledge on the phase structure is however still lim- 
ited and the description of the matter around the phase 
transitions does not reach a consensus, where a typical 
size of the critical temperature and chemical potential is 
considered to be of order Aqcd- The phases of QCD are 
characterized by symmetries and their breaking pattern: 
QCD at asymptotically high density leads to the color- 
flavor-locked phase as the true ground state under the 
symmetry breaking pattern, SU{3)c x SU{3)l x SU{3)ii 
down to the diagonal subgroup SU{3)c+l+r. [2]. The 
residual discrete symmetries characterize the spectra of 
excitations. 

At zero temperature and density, an alternative pat- 
tern of spontaneous chiral symmetry breaking was sug- 
gested in the context of QCD [1, |3, H. This pat- 
tern keeps the center of chiral group unbroken, i.e. 
SU{Nf)L X SU{Nf)R SUiNf)v x (ZatJa, where a 
discrete symmetry {Z]\Jj:)a is the maximal axial subgroup 
of SU{Nf)L X SU{Nf)R. The symmetry protects a 
theory from condensate of quark bilinears (qq). Spon- 
taneous symmetry breaking is driven by quartic conden- 
sates which are invariant under both SU{Nf)v and Zj^j, 
transformation. Although meson phenomenology with 
this breaking pattern seems to explain the reality reason- 
ably [3j] , this possibility is strictly ruled out in QCD both 
at zero and finite temperatures but at zero density since 
a different way of coupling of Nambu-Goldstone bosons 
to pseudo-scalar density violates QCD inequalities for 
density-density correlators |6j. However, this does not 
exclude the unorthodox pattern in the presence of dense 
baryonic matter. There are several attempts which dy- 
namically generate a similar breaking pattern in an 0(2) 
scalar model [?] and in A/" = 1 Super Yang-Mills the- 



ory [S|J. 

Within the Skyrme model on crystal, a new intermedi- 
ate phase where a skyrmion turns into two half skyrmions 
was numerically found This phase is character- 

ized by a vanishing quark condensate (qq) and a non- 
vanishing pion decay constant. Recently, another novel 
view of dense matter, Quarkyonic Phase, has been pro- 
posed based on the argument using the large Nc counting 
where Nc denotes the number of colors |10| : In the large 
Nc limit there are three phases which are rigorously dis- 
tinguished using the Polyakov loop expectation value ($) 
and the baryon number density (Nb)- The quarkyonic 
phase is characterized by (4>) = indicating the system 
confined and non- vanishing (Nb) above fJ-B — Mb with a 
baryon mass Mb ■ The separation of the quarkyonic from 
hadronic phase is not clear any more in a system with fi- 
nite Nc- Nevertheless, an abrupt change in the baryon 
number density would be interpreted as the quarkyonic 
transition which separates meson dominant from baryon 
dominant regions. This might appear near the bound- 
ary for chemical equilibrium at which one would expect a 
rapid change in the number of degrees of freedom 

A steep increase in the baryon number density and 
the corresponding maximum in its susceptibility xb are 
driven by a phase transition from chirally broken to re- 
stored phase in most model-approaches. Interplay be- 
tween (de) confinement and chiral symmetry breaking 
has been studied within a Nambu-Jona-Lassinio model 
with Polyakov loops which describes how the decon- 
finement and chiral phase boundaries are changed from 
Nc = 00 down to A^c = 3 Ts*!. The model study shows 
that the chiral phase transition at T = appears just 
above mass threshold hb = Mb and thus a large xb is 
associated with the chiral phase transition. However, a 
constituent-quark picture does not directly describe the 
thermodynamics of hadronic matter and there are no 
a priori reasons that the quarkyonic transition should 
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be accompanied by chiral phase transition. Besides, it 
seems unlikely that the chiral symmetry is (even par- 
tially) restored slightly above the freeze-out curve where 
the baryon density is not high enough to drive a phase 
transition. From this perspective, further investigations 
of dense baryonic matter and a possible appearance of the 
quarkyonic phase in QCD with Nc = 3 require a model- 
ing in terms of dynamical hadronic-degrees of freedom in 
a systematic way. 

In this paper we will address this issue under the al- 
ternative pattern of chiral symmetry breaking in dense 
hadronic matter. We will show a possible intermediate 
phase between chiral symmetry broken and its restored 
phases with analyses using a general Ginzburg-Landau 
free energy. This leads to multiple critical points and one 
of them is associated with restoration of the center sym- 
metry rather than that of chiral symmetry. In the new 
phase with unbroken center symmetry the net baryon 
number susceptibility exhibits a strong enhancement al- 
though the chiral symmetry remains spontaneously bro- 
ken. This is reminiscent of the quarkyonic transition and 
our framework provides a theoretical description of the 
quarkyonic phase on the bases of a chiral Lagrangian with 
two distinct order parameters. An analysis using a linear 
sigma model for hadron mass spectra is also made. 



2. A MODEL FOR 2-QUARK AND 4-QUARK 
STATES 



We construct a chiral Lagrangian for 2- and 4-quark 
states under the following pattern of symmetry breaking. 



S /2. The M and E are expressed as 



where a and x represent scalar fields and <p and ip pseu- 
doscalar fields, and e^fe is the total anti-symmetric tensor 
with £123 = 1. In general the field S contains an isospin 
2 state. One can take appropriate parameters in a La- 
grangian in such a way that this exotic particle is very 
heavy. Thus, we will consider only isospin (x) and 
1 (ip) states in this paper. The fields transform under 
SU{2)l X SU{2)ii as chiral non-singlet. 



M 



„(3)y.„(3)t 



(2.4) 



This transformation property implies that the field M 
changes its sign under the center Z2 of SU(2)i (or 
SU(2)fl), while S is invariant: 



M ■ 



-M , 



(2.5) 



Up to the fourth order in fields one obtams a potential, 

2 \ 2 

V{M,J:) = -^Tr [MM^] + y (Tr [MM^]f 

\ T. \ 1 T ^ 1 V 1 ^ ^2"" . , ^ /j-i "^2 / V 1 ^ ^ \ 2 

+ 2giEahTr [TaMTtM^ + gzEafcS^.Tr [MM^ 

+ g3DetT. + gi(DetM + h.c.) . (2.6) 



SU{Nf)L X SU{Nf)i 



SUiNf)v X {ZNf)A 

SU{Nf)v. (2.1) 



In this paper we will restrict ourselves to a two-flavor 
case. 



The last term violates the U{1)a symmetry. The coeffi- 
cients of the quartic terms are positive for this potential 
to be bounded. Other parameters gi can be both posi- 
tive and negative and will determine the topology of the 
phase structure. An explicit chiral symmetry breaking 
can be introduced through, e.g., 



2.1. Lagrangian 

We introduce a 2-quark state M in the fundamental 
and a 4-quark state T, in the adjoint representation as 



Mij ~ qR,jqL,i , 



Sab 



qLTalt,qLqRTbT qR ■ 



(2.2) 



Vsb{M,^) = -h<j-ah\, 



(2.7) 



with constants h and a. Note that a similar Lagrangian 
was considered for a system with 2- and 4-quark states 
under the symmetry breaking pattern without unbroken 
center symmetry in [l^ where their 4-quark states are 
chiral singlet and the potential does not include quartic 
terms in fields. 



where the flavor indices run = 1,2 and {a,b,c) = 

1,2,3 and Pauh matrices = 2r° with tr[T'^T''] = 



2.2. Ginzburg-Landau effective potential 

We first study possible phases derived from the effec- 
tive potential (|2.6p taking 



We consider E as any linear combination of qq-qq and qq-qq type 
fields allowed by symmetries. 



Mij = -j=(j5ij , Eab 
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FIG. 1: Phase diagram with D = F = and /i = 0. 
The sohd and dashed lines indicate first and second or- 
der phase boundaries, respectively. One tricritical point, 
TCPi, is located at {A,B) = (0,1/4) and another, TCP2, 
at {A, B) = (1/4, -1/8). The triple point represented by T is 
at \a,B) = (1/8,0). 



One can reduce Eq. (|2.6p as well as an explicit breaking 
term to 



(2.9) 



We will take C = — 1 without loss of generality in the 
following calculations. 

We start with the potential for _D = F = and ft, = 0, 



V = Aa^ + Bx" 



X 



(2.10) 



Phases from this potential can be classified by the coeffi- 
cients A and B. The expression of the phase boundaries 
is summarized in Appendix \^ Here we discuss the ob- 
tained phase structure shown in Fig. [T] There are three 
distinct phases characterized by two order parameters: 
Phase I represents the system where both chiral sym- 
metry and its center are spontaneously broken due to 
non- vanishing expectation values xo and co- The center 
symmetry is restored when ctq becomes zero. However, 
chiral symmetry remains broken as long as one has non- 
vanishing indicated by phase II. The chiral symmetry 
restoration takes place under Xo which corresponds 
to phase III. The phases II and HI are separated by a 
second-order line, while the broken phase I from II or 
from III is by both first- and second-order lines. Accord- 
ingly, there exist two tricritical points (TCPs) and one 
triple point. One of these TCP, TCP2 in Fig. [H is as- 
sociated with the center Z2 symmetry restoration rather 
than the chiral transition. 

Two phase transitions are characterized by suscepti- 
bilities of the corresponding order parameters. We intro- 
duce a 2-by-2 matrix composed of the second derivatives 
of V as 



C 



(2.11) 



with 



dx^ 



92 y 

r — c — 



(2.12) 



under the solutions of the gap equations, ctq and xo- A 
set of susceptibilities is defined by the inverse of C [IQ] ; 



1 



X = 



detC 



(2.13) 



We identify the susceptibilities associated with 2-quark 
and 4-quark states as 



X2Q — Xii 7 X4Q — X22 



(2.14) 



The X2Q is responsible to the Z2 symmetry and the XiQ 
to the chiral symmetry restoration. 

We consider X2Q and X4Q around the TCPi in Fig. [T] 
where the potential has zero curvature and thus detC = 
0. When approaching the TCPi from broken phase I by 
tuning A and B e.s A ^ ^critical ~ and B — 1/4, these 
susceptibilities diverge as 



X2Q 



X4Q 



t 



-2/3 



(2.15) 



where Acriticai — A ^ t with the reduced temperature 
or chemical potential, e.g. t = \fi — ^c|/Mc- The gap 
equations determine the scaling of 2-quark and 4-quark 
condensates as 



1/3 



Xo ^ t 



1/3 



(2.16) 



Consequently, the quark number susceptibility Xq 
—d^V/dfi^ exhibits a singularity as 



Xq 



■ X2Q t 



-2/3 



(2.17) 



This critical exponent is same as the one in the 3-d Ising 
model. The coincidence can be understood due to the 
same Z2 symmetries 

The critical behavior near the TCP2 involves more: 
When the A is approached as 1/4 — t with B = —1/8 
fixed, X2Q and X4Q diverge as 



X2Q 



X4Q 



t 



-1/2 



(2.18) 



and only ao vanishes as CTq ^ t^/"^. As a result, the quark 
number susceptibility Xq diverges as 



Xq 



-1/2 



(2.19) 



The Z2 symmetry in the 3-d Ising system is not the center of 
two-flavor chiral group, but emerges in the direction of a linear 
combination of quark number and scalar densities (16|| . 
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Note that the critical exponent 1/2 is different from the 
one near the TCPi, which may reflect different symme- 
tries possessed by the system at TCP2, SU{2)v and the 
center Z2, from that at TCPi, SU{2)l x SU{2)r includ- 
ing its center {Z2)l x (^2)7?,- Those exponents at TCPi, 2 
are changed when D ^ (see below). 

When the second-order phase transition separating 
phase I from 11 or from 111 is approached from the broken 
phase with a fixed B, we have 



X2Q 



X4Q 



1 

B 



(2.20) 



where B is a finite number, which thus gives no singular- 
ities in X4Q- The 2-quark condensate scales as CTq 
and the quark number susceptibility Xq is finite along the 
second-order phase transition line: 



Xq 



(2.21) 



Nevertheless, Xq is enhanced toward the phase transi- 
tion induced by X2Q and becomes small above the phase 
transition. Such abrupt changes in Xq indicate the phase 
transition, especially for a negative B which is driven by 
the center symmetry restoration rather than the chiral 
phase transition. 

Near the second-order chiral transition between phase 
II and III, one obtains from B ^ t 



Xo 



X4Q 



(2.22) 



Since the chiral symmetry including the center symmetry 
prohibits the Yukawa- type coupling of x to a fermion and 
an anti-fcrmion in the fundamental representation, the 
coupling of X to the baryon number current would be 
highly suppressed. Therefore, Xq shows less sensitivity 
around the chiral transition. 

Once small h is turned on, chiral symmetry and its cen- 
ter are explicitly broken. Second-order phase boundaries 
are replaced with cross over and the two TCPs with two 
critical points. The singularity in Xq is now governed by 
the Z2 universality class of 3-d Ising systems. Thus, the 
scaling of Xq at the critical points (CPs) will be given by 



Xq 



t 



-2/3 



(2.23) 



A cubic term in x modifies the previous phase structure 
shown in Fig. [TJ The phase diagram from the potential. 



V ^Aa^ + Bx^ 



a^ + x' 



Dx^ 



(2.24) 



is classified by the following regions of D: (i) — 1 < I? < 
0, (ii) D < -I, (iii) < D < 1 and (iv) I < D. One 
observes a deformation of the boundary lines depending 



'^■^ As we will show below, the phase transition from phase II to 
phase III is of first order in a more general parameter choice. 
Thus, Xq exhibits a jump at the chiral phase transition point. 





CPi 


TCPi 


TCP2 


D< 


2/3 




1/2 


D = 




2/3 


1/2 


D > 




1/2 





TABLE I: The critical exponents of the quark number sus- 
ceptibility for vanishing and non-vanishing D at two tricritical 
points and at the critical point (CP). 



on D as in Fig. [51 The phase transition line separating 
phase II from phase III becomes of first order due to the 
presence of . For any negative D, (i) and (ii), a crit- 
ical point CPi appears as a remnant of TCPi for D — Q. 
TCP2 remains on the phase diagram for — 1 < D < 0, 
(i), which eventually coincides with the triple point at 
D = —1, (ii). For positive D, (in) and (iv), the transi- 
tion line which separates phase 1 from phase II turns to 
be of first order everywhere. The triple point approaches 
the TCPi and coincides when a positive D reaches unity. 
The different order of phase transition between phase 1 
and phase II for -1 < Z) < to that for < D < 1 
can be understood as follows: For D = Q (see FigH]) the 
vacuum expectation value (VEV) xo is positive in phase 
I near the phase boundary between phase I and II due to 
the existence of the — ct^x term in the potential. In phase 
II, on the other hand, when the positive xo provides a 
local minimum of the potential, — xo also does, and both 
coincide with the global minima. These two vacua are 
physically equivalent, so that the phase transition from 
phase I to phase II can be of second order. When we add 
Dx^ term with negative D to the potential, the local min- 
imum corresponding to the positive xo is only the global 
minimum in phase 11. This can be smoothly connected to 
the vacuum in phase 1 where the VEV xo is positive. On 
the other hand, when D is positive, the negative xo gives 
the global minimum in phase II. Thus, there is a mis- 
match of Xo along the phase boundary separating phase 
I from phase II, which indicates a first-order transition. 

D also affects the quark number susceptibility Xq- ^s 
in the case of -D = 0, the Xg exhibits a more relevant 
increase toward the Z2 symmetry restoration than at the 
chiral phase transition. The critical exponents of Xq is 
summarized in Table HI One finds that the two regions, 
D < Q and < Z?, corresponds to different universality. 
The cubic term plays a similar role to an explicit sym- 
metry breaking term in the potential. This may be an 
origin for the appearance of a critical point. 

For — 1 < D < 0, TCP2 for h = becomes a critical 
point, CP2, for finite h. When the value of h is increased, 
the CP2 approaches the triple point and coincides with 
it for a certain value of ft., /iq. The topology of the phase 
diagram for larger h > hg agrees with that for D < — 1. 
Similarly, the TCPi in the < D < 1 phase diagram 
becomes a critical point CPi and disappears for a suffi- 
ciently large h. On the other hand, the CPi stays in the 
phase diagram Fig. [2](i) and (ii) for any value of h. The 
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FIG. 2: Phase diagram for different values of D under F = and h = Q. The sohd and dashed hnes indicate first and second 
order phase boundaries, respectively. 



scaling of Xq there will be given by 

Xq - i"'/' • (2.25) 

We note that adding finite F to the potential does not 
generate any essential differences from the above result 
with F = 0. 



3. HYPOTHETICAL PHASE DIAGRAM AND 
QUARK NUMBER SUSCEPTIBILITY 

From the above observations one would expect phase 
diagrams mapped onto (T, fi) plane, fn the chiral limit 
a new phase where the center symmetry is unbroken but 
chiral symmetry remains broken might appear in dense 
matter since a.t fi — this phase is strictly forbidden by 
the no-go theorem. With an explicit breaking of chiral 
symmetry one would draw a phase diagram as in Fig. [31 
The intermediate phase remains characterized by a small 
condensation |(To| ^ |xo|- One would expect a new crit- 
ical point associated with the restoration of the center 
symmetry, CP2, rather than that of the chiral symmetry 
if dynamics prefers a negative coefficient of the cubit term 



in X- Multiple critical points in principle can be observed 
as singularities of the quark number susceptibility. 

It has been suggested that a similar critical point in 
lower temperature could appear in the QCD phase di- 
agram based on the two-flavored Nambu-Jona-Lasinio 
model with vector interaction [17] and a Ginzb urg - 
Landau potential with the effect of axial anomaly [18| . 
There the interplay between the chiral (2-quark) con- 
densate and BCS pairings plays an important role. In 
our framework without diquarks, the critical point dis- 
cussed in Fig. [3] (left) is driven by the interplay between 
the 2-quark and 4-quark condensates, and is associated 
with restoration of the center symmetry where anomalies 
have nothing to do with its appearance. Nevertheless, 
the cross over in low temperatures may have a close con- 
nection to the quark-hadron continuity [IS] and it is an 
interesting issue to explore a possibility of dynamical cen- 
ter symmetry breaking in microscopic calculations. The 
present potential (|2.9p leads to a first-oder transition of 
chiral symmetry even with an explicit breaking. This 
may be replaced with a cross over when one considers 
higher order terms in fields and other symmetry break- 
ing terms as well as in-medium correlations to baryonic 
excitations, which is beyond the scope of this paper. 
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FIG. 3: Schematic phase diagram mapped onto (T, fi) plane with a negative D (left) and with a positive D (right). The solid 
lines indicate first order phase boundaries, and dashed lines correspond to cross over. 
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FIG. 4; The behavior of the baryon number susceptibility as 
a function of chemical potential assuming the phase diagram 
of Fig. [3] (left). The condensates and the susceptibility show 
a jump also at /iz2 when the phase structure of Fig. [3] (right) 
is preferred. 



FIG. 5: A schematic behavior of the susceptibility X2Q near 
the CPi as a function of temperature assuming the phase 
diagram of Fig. [3] (left) . 



Appearance of the above intermediate phase seems to 
have a similarity to the notion of Quarkyonic Phase [lol . 
[l3| , which is originally proposed as a phase of dense mat- 
ter in large Nc limit. The transition from hadronic to 
quarkyonic world can be characterized by a rapid change 
in the net baryon number density. This feature is driven 
by the restoration of center symmetry and is due to the 
fact that the Yukawa coupling of x to baryons is not al- 
lowed by the Z2 invariance. Fig.[l]shows an expected be- 
havior of the quark (baryon) number susceptibility which 
exhibits a maximum when across the Z2 cross over. This 
can be interpreted as the realization of the quarkyonic 
transition in Nc — 3 world. How far fi22 from /ichirai is 
depends crucially on its dynamical-model description. 



Thus, the present analysis does not exclude the possibility that 
both transitions take place simultaneously and in such case en- 
hancement of xb is driven by chiral phase transition. The phase 
with Xo 7^ and ao = does not seem to appear in the large 



It should be noticed that the critical point in low den- 
sity region, indicated by CPi in Fig. [3] (left), is different 
from a usually considered CP [2^ in the sense that the 
CPi is not on the cross over line attached to the T = 
axis. When we take a path from the broken phase (both 
(To and Xo are large) to the symmetric phase (both ctq and 
Xo are small) passing near the CPi, the X2Q may exhibit 
two peaks; one is located near CPi and another is on the 
cross over line. We show a schematic behavior of X2Q as 
a function of temperature, together with cto and xo hi 
Fig. [5l The appearance of two peaks in X2Q reflects the 
fact that (To becomes small across the CPi and the cross 
over. The first decrease in (To near CPi is caused by a 
dropping xo, while the second is by the chiral symmetry 
restoration. 



Nc Umit [1, 0, 0I- It would be expected that including 1/Nc 
corrections induce a phase with unbroken center symmetry. 
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4. HADRON MASS SPECTRA AND PION 
DECAY CONSTANT 

In this section we derive meson mass spectra in a linear 
sigma model. The Lagrangian with the potential (|2.6p is 
expressed in terms of the mesonic fields as 



Obviously, the determinant of the above mass matrix for 
(j) and "0 is zero and thus massless pseudo-scalar fields are 
a mixture of 2-quark and 4-quark states. 

The mass eigenstates are introduced with a rotation 
matrix as 



cos a sm 
- sin 6 cos 



with 

U - 




cos f sm 
— sin 9 cos 




(4.8) 



m 
"2 



2 

A2 



with the angles 



-'^ 4 I 2 2 7*2 

3X +-^x^ 



1 / ^ N 2 



tan (26*) 



2V^ 



9 9 

— 



tan(20)^i^5^. (4.9) 



V3 U ^2^^ 



The masses of scalar mesons are give by 



1 



V2a(j)-'if 
(4.2) 



TO5 = TO^ cos^ 6* + sin^ 6* — Vs^ctq sin(26') , 

= rn2 j,Qg2 5) _,_ gjj^2 5/ + ^gcTQ sin(26') , (4.10) 



where gi = —g {g > 0) and 52 = were taken. In 
addition, we set 54 = since the (74-term generates only 
a shift in for A^/ = 2. We also set the explicit breaking 
being zero. 

The condensate of the mesonic fields in the phase 
where both the chiral symmetry and its center Z2 are 
broken are determined from the coupled gap equations 
given by 



and those of pseudo-scalar mesons by 

.9(3^2 + 8x1) 



mp = , Trip, ~ 



with 



\/3ag + 8x1 



sin 9 



2^3X0 ' 



(4.11) 



(4.12) 



V3.g 
1 



^2 -2 , 53 
y Xo -m + -^xo ] Xo : 



The pion decay constant is read from the Noether cur- 
rent, - aodt'cj) + 4/^6x09^-0, as 



Xo 



(4.3) 
(4.4) 



V3g 

with = A? + 3A|. Shifting the fields as 
cr->a-|-cro, X^X + Xo, 
the potential reads 

1 22|1 2 7'2|1 22,1 2/2 

— \figaQ ax ~ V^ga^ • ^ + • • • , (4.5) 



gXg. 



(4.13) 



U 



where ellipses stand for the terms including the fields 
more than three, and 



2 o \ 2 2 2 

m„ = 2A CTo , 711 = 



9 2 2;-2 2 , 



2 xo 
V3 



1 

3' 



2 Xo 



(4.6) 



Since we consider a system in the chiral limit, the massive 
P' state is decoupled from the current and F.^' =0, as it 
should be. It should be noted that, when |cro| ^ |xo|, the 
NG boson is dominantly the 2-quark state. The 4-quark 
component becomes more relevant for \/3|(To| < V8|xo|j 
i.e. 9 > 7r/4. 

When the coupling 173 is negative, which corresponds 
to D < in the Ginzburg-Landau potential given in sec- 
tion [21 the phase transition from phase I (ctq ^ and 
Xo 7^ 0) to phase II (ctq = and xo 0) can be of second 
order. In such a case, the restoration of the center Z2 
symmetry is characterized by vanishing ctq ■ Approaching 
the restoration from broken phase, one finds the lowest 
scalar meson mass degenerate with the P state, while the 
pion decay constant remains finite due to Xo 7^ 0; 



The mass terms thus become 

n2 



ms ^ mp = , Ft, 



■Xo, 



2 \-V3gao mf. 



1 



ml -V2g(To 



-V2 



50-0 




with 



(4.7) 



Xo 



3to2 / \/S 



93 



2A2 



V3 



2A2 



53 



(4.14) 



(4.15) 
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The vanishing S'-state mass corresponds to a diver- 
gence of the susceptibiUty X2Q, which is responsible to 
restoration of the center symmetry. The scalar S and 
pseudo-scalar P states thus become the chiral partners 
on the phase boundary. In the Z2 symmetric phase the 
meson masses are found from the potential (14. 2|) as 

ml ^ -m'^ - gXo , = -rn^ + -^Xo , 

ml^^~X\l + ^Xo, m| = 0. (4.16) 

There is no mixing in this phase, tan0 — tanO — 0, so 
that fj, (j>, Xi '4' ^re the mass eigenstates. This implies 
that the pure 4-quark state is the massless NG boson 
in the Z2 symmetric phase. Due to the broken chiral 
symmetry, a and (f> states are not degenerate in mass. 
The vector and axial-vector states neither degenerate in 
mass @, since both vector and axial- vector currents are 
invariant under the Z2 transformation but broken chiral 
symmetry does not dictate the same masses. 

When I (73/5 1 ^ 1, the chiral phase transition from 
phase II ((To = and xo 7^ 0) to phase III (ctq = and 
Xo = 0) will be of weak first-order. In this case, xo then 
Ft, approach zero near the phase transition point. This 
is controlled by m approaching zero, which corresponds 
to B approaching zero in the Ginzburg-Landau poten- 
tial discussed in section 12.21 The iso-spin 2 state will 
become very light near the phase transition. This may 
suggest that, when the gsDetE term is small and the chi- 
ral phase transition is of weak first-order, a light exotic 
states with 1 = 2 might exist in dense baryonic matter. 
When there exists the non-negligible g^DefE term, on 
the other hand, such state never becomes light since the 
chiral phase transition is of strong first-order. 

In two flavors, the system would prefer the parity dou- 
bling for baryons in the Z2 symmetric phase where the 
VEV Xo does not generate the baryon masses [(|. In 
the parity doubling scenario [2l|, all the baryons have 
their parity partners and then each pair of parity part- 
ners has a degenerate mass. On the other hand, in the 
naive scenario the lightest baryon does not have a parity 
partner, so that it becomes massless in the Z2 symmetric 
phase. We list hadron mass spectra expected in phase I 
and phase II in Table HIl 



When we approach the phase boundary from the Z2 symmetric 
phase to the Z2 broken phase, in Eq. 114.161 1 approaches zero, 
since —m? = VSffXO is satisfied at the phase boundary. The 
pseudoscalar mass approaches -^gxo which coincides with 
the mass of P' in the Z2 broken phase (see Eq. | |4.11| |). 

In reference Q the degeneracy of the massive scalar and pseu- 
doscalar mesons made of 4-quarks carrying the same isospin for 
a general number of flavors was shown. In case oi Nf =2 the 
U(1)a anomaly generates a mass difference between the a state 
and the pseudoscalar meson with / = {rj). In the present anal- 
ysis, we did not include the 7 = pseudoscalar and the 7 = 1 
scalar mesons from the beginning by assuming that they are very 
heavy. 



5. CONCLUSIONS 

We have discussed a new phase where chiral symme- 
try is spontaneously broken while its center symmetry is 
restored. This might appear as an intermediate state be- 
tween chirally broken and restored phases in (T, fi) plane. 
The appearance of the intermediate phase with unbroken 
Z2 also suggests a new critical point associated with the 
center symmetry in low temperatures. A tendency of the 
center symmetry restoration is carried by the net baryon 
number density which shows a rapid increase and this is 
reminiscent of the quarkyonic transition. In this phase 
there exist the NG bosons and thus nothing prevents the 
anomaly matching. It has been shown that the anomaly 
matching conditions are also valid in gauge theories at 
finite density [13] . This may suggest that the phase III is 
chirally restored and deconfined. The U(1)a symmetry 
remains broken and the heavy 77 mass can be controlled 
with a certain anomaly coefficient. 

There are subtleties in baryon masses since the exis- 
tence of the center symmetry does not immediately dic- 
tate the parity doubling for a general number of fiavors: 
Here we consider the case in massless three fiavors. The 

53- term in p.6p now generates x* contribution, while the 

54- term does (T^ one. It follows that Dx"^ is removed from 
(|2.9p and another cubic term is added. Omitting the 
cubic term results in the same phase diagram as Fig.[T] 
with two TCPs. When the cubic term is included, it 
is conceivable that phase II and phase III in Fig[2] are 
separated by a second-order phase boundary, which will 
become a first-order one when we take quantum fluctua- 
tions into account ^23*1 . The topologies are expected to be 
quite similar to those shown in Fig. [2l so that we expect 
a strong enhancement of the quark number susceptibility 
at the Z3 restoration point. Differently from the case for 
Nf = 2, the Y.ab field is allowed to couple to the octet 
baryon states as, e.g. Ba^abBb, and the baryon number 
current couples to the x state which becomes massless at 
the chiral restoration point. As a result, the quark num- 
ber susceptibility might show another peak at the chiral 
restoration. Hadron masses in the Z3 symmetric phase 
are slightly different from those under Z2 invariance: In 
the mesonic sector the party partners are degenerate and 
the degeneracy does not generally occur in the baryonic 
sector 0. Following Ref. Q, possible operators for the 
baryons are expressed as 

Bi = {qLqLqR)L , Bl = {qLqRqR)L , BI = {qLqLqL)L , 
Br = {qRqRqL)R , B\ = {qRqLqL)R , = {qRqRqR)R , 

(5.1) 

where the color and flavor indices are omitted. For 
the octet baryons, the representations under the chiral 
SU(3)LxSU(3)fl of these baryonic fields are assigned as 

i?i^(3,3), Si ^(3, 3), Bi^(8,l), 
S]j^(3,3), i?2^(3,3), i?i^(l,8). (5.2) 

When the is the lightest octet baryon, which we call 
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phase I: cro / , xo / 


phase II; (70 = , xo / 


SU{2)v 


SU{2)v X (Z2)a 


ms / , mp — 


ms7^mp/0,mp/ = 


mv / mA 


mv / rriA 


F. = ^<tI + (8/3)xg 






1 mjY+ = O(ground state) 


mjv+ / 


(i) naive: < m^,+ = m^,/- ^ 
[ (excited states) 




/■•\ • J "T-]v+ = "ijv- / 
(u) mirror: s 

1 (all states) 





TABLE II: The mass spectra of mesons and baryons in different phases for Nf = 2. Baryons transform with the naive chirality 
assignment as V'H.i ~^ gR,L^R,L, while with the mirror assignment as ipiR.L —> gR,Li^iR.L and '4'2r,l gL.Ri'2R,L with 
gR,L G SU{2)r^l where two nucleons tpi and ip2 belong to the same chiral multiplets. 



the naive assignment, it is still massive in the Z3 sym- 
metric phase, since the Yukavi^a coupling of the 4-quark 
state Tiab is possible as, e.g. BaX'abBb- When the Ughtest 
baryons are described by a combination of and i?^, 
which we call the mirror assignment, they are degenerate 
with each other in the Z3 symmetric phase. We summa- 
rize these features in Table IIIII The baryon masses cru- 
cially depend on a way of chirality assignment. It would 
be an interesting issue to clarify this within a more elab- 
orated model. 

The main assumption in this paper is a dynamical 
breaking of chiral symmetry SU{Nf)L x SU{Nf)ii down 
to a non-standard SU{Nf)v x (Zj^^)^ although this 
seems to be theoretically self-consistent. Calculations us- 
ing the Swinger-Dyson equations or Nambu-Jona-Lasinio 
type models with careful treatment of the quartic op- 
erators may directly evaluate this reliability. Besides, 
anomalously hght NG bosons, ~ O(m^), could lead 
to an s-wave pion condensation as discussed in A 
calculation using the Skyrme model shows a similar in- 
termediate phase 0. Although the above non-standard 
pattern of symmetry breaking was not imposed in the 
Skyrme Lagrangian, the result could suggest an emergent 



J 



symmetry in dense medium. This intermediate phase 
would be an intriguing candidate of the quarkyonic phase 
if it could sustain in actual QCD at finite density and 
would lead to a new landscape of dense baryonic matter. 
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APPENDIX A: PHASE BOUNDARIES FROM GINZBURG-LANDAU POTENTIAL 



The relevant expressions for the phase boundaries obtained from the potential (|2.9p are given below. We will take 
D = F = and the chiral limit h = 0. 



• Second-order phase transition when B > 1/4: 

The solutions for a and x on this boundary are given by 

(ao,xo) = (0,0). 

• First-order phase transition when < B < 1/4 and < A < 1/8: 



A = 




(A.l) 
(A.2) 

(A.3) 
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phase I: cro / , xo / 


phase II: ao — ,xo ^ 


SU{Nf)v 


SU{Nf)v X {ZNf)A 


ms / , mp — 
mv / niA 


ms = mp / 0, mp/ = 
mv 7^ JTiA 


mjv+ / 


(i) naive: mjv+ 7^ 

(ii) mirror: mj^+ = m^y- 7^ 



TABLE III: Same as in Table |ll] but for Nf = 3. 



The solutions for a and x on this boundary are given by 
('TcXo) = (0,0), 











± 













1 f 1 



3 / 1 



5\ 




• First-order phase transition when — l/8<i3<0 and 1/8<j4<1/4: 

.-.is. 

The solutions for a and x on this boundary are given by 



Second-order phase transition when B < —1/8 and A> 1/4: 



The solutions for a and x on this boundary are given by 



(cto ,Xo) = 0,± 



-B 



• Second-order phase transition when A > 1/8: 

B = 0. 

The solutions for ct and x on this boundary are given by 

(ao,Xo) = (0,0). 

I 



(A.4) 



8 2\ 4, 



1/2 



(A.5) 



(A.6) 



(A.7) 



(A.8) 

(A.9) 

(A.IO) 
(A.ll) 
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